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ABSTRACT 


The objective of this thesis was to develop a computer 
program for use in evaluating data acquired irom the)Weparv= 
ment of Aeronautics Cessna 310H model aircraft. The program 
was designed for use with the IBM 360 computer, in conjunc-— 
tion with a CALCOMP plotter. The program format is very 
similar to that of IBM’s Scientific Subroutine Package 
routines; in addition, the program provides the student with 
a detailed set of instructions for loading the data. 

Using information taken from the engine curves for the 
aircraft, the program evaluates input data, determines the 
power required for flight at altitude, reduces the values 
to equivalent sea level values and outputs a graph of power 
versus velocity for level flight. Lift and drag coefficients 
are calculated as functions of aircraft weight, velocity, 


and power required, and are plotted on a second graph. 
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LT. INTRODUCTION 


The laboratory course is one of the most demanding for 
the graduate student because of the extremely time consuming 
process of reducing raw data. Often the student may find 
himself neglecting important course work because of the 
large amount of work required to reduce laboratory data by 
hand to meaningful results. This is, of course, an unsat- 
isfactory situation, but is one which seems to be the rule 
rather than the exception for many students. 

The Department of Aeronautics Cessna 310H aircraft 
flying laboratory, which supports the Flight Evaluation 
Techniques courses, AE 3321 and AE 3322, at the Naval Post- 
graduate school, is typical of such time consuming labora- 
tories. With the incorporation of a photo-panel data 
acquisition system in the aircraft came the capability to 
collect large amounts of in-flight data which require large 
amounts of the student’s time for evaluation. 

There appeared to be a need for some sort of procedure 
which would reduce the student workload in this laboratory 
by relieving him of the burden of hand reducing the large 
volume of raw data and hand plotting the results. A program 
which would require Sale that the student assemble the data 
in a proper sequence and input it to a computer, which in 
turn would evaluate the data and SubpuD tabular and sraphical 
displays of the results, would be extremely beneficial in 


lightening the student workload. Such a program would free 


ee = , “ee 





the student, enabling him to spend time evaluating his 
results instead of tabulating them. With this need in mind 
aecompuLer propram tor the reduction of in-flight data for 
the Cessna 310H aircraft was written. 

CESSNA was designed as a student aid, with the intent 
that it be self-explanatory, containing all information and 
Uistructacens Necessary for the studentruio run the program. 
Parts II and III of this thesis may be used as a user’s 
manual in the event more detailed information is required, 
Sumitenodit edt oneand, oreaddiciens arerdiesired. Thesprogzram 
is written in FORTRAN language but knowledge of the language 
is not required, other than the ability to type data cards 
and a job card in proper format; this format information is 
readilv available at the computer center. The card deck is 
complete. The student needs only to add his data cards to 


the end of the deck and it is ready to run. 





IIT. PROGRAM ORGANIZATION AND OPERATION 


A. GENERAL ORGANIZATION 
As shown in Fig. 1, the complete program consists of 


two primary sections and two subroutines: 


iPeeiccrare Power fequirementes. 

£e Lift and Drag: 

3. Interpolation Subroutine: 

4. Sorting Subroutine. 

When program execution begins, aircraft engine data are 
read in for internal use in the program. Following this, 
student data concerning aircraft weight, altitude, velocity, 
RPM and manifold pressure, and temperature, are read in for 
evaluation. 

CESSNA uses the engine data and the student’s RPM and 
manifold pressure data to determine shaft horsepower required 
for level flight at the operating altitude for the input 
velocity range. Horsepower, Bitoni with velocity and weight, 
is then used to determine lift and drag coefficients for 
the aircraft. Results are printed out in tabular form and 


are eraphed using a CALCOMP plotter. 


Bs ATRCRAFT POWER REQUIREMENTS 

The heart of the program is the set of engine perfor- 
mance curves shown in Fig. 2. On the sea level side of the 
chart the curves are a linearized version of the original 


curves (Continental Motors Corporation, 1958), and are a 


a4 Y 
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Figure 1 


General Program Flow Chart 
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good approximation to the originals. In the recommended 
cruising range the maximum error in the linearized curves 
is approximately 1 to 1.5% in the mid-range RPM’s and mani- 
fold pressures (or a maximum of approximately 2.5 brake 
horsepower), and much lower in the upper and lower ranges, 
where the original curves are very nearly linear. The 
altitude curves are a copy of the origina& curves. 

The curves are utilized in the following manner. For 
pour hth lcmic Smee Scay level scleolmeurves, Uhegslope of 
the line, maximum manifold pressure, and maximum brake horse- 
power have been determined. For the altitude set of curves 
each RPM line has had slope and maximum brake horsepower 
determined. A question may arise at this point as to why, 
for a given RPM, the maximum brake horsepower is not the 
same for both curves and why, therefore, one horsepower 
valve Mer noU=suiticient for any calculations needed. The 
answer is two-fold. First, in linearizing the sea level 
set of curves, the maximum points for the curves changed, 
whereas the maximum points for the altitude side did not; 
hence, two values. Secondly, since points on the RPM lines 
are determined using slopes and maximum points, it is 
necessary to keep both values so that calculations for both 
sets of curves will be accurate. | 

A simplifying assumption was made concerning the mani- 
fold pressure lines on the altitude side. For purposes of 
this program the lines were assumed to be vertical. This 
introduces an error of not more than 1.1%, which yields at 


most a.difference of two brake horsepower. This assumption 


Ze 

















€l-1 TAO 
160) 8 OOO! JON LILA 
Oz SI Ol 68GH2 9 Sieg Zz |) oO 
U 
(D) OILVY ALISN3A 
SO 90 LO 80 6'0 O'1 





dvb! JONVANNOSH3d Baniilay O09 


Te) 
~O 
i 
re) 
NA 
a 
N 
NI 
Oo Oo 
oO re) 


mm 

N 
4 
N 


C/O 
ie 
uw) 
WN 
© 
WN 
a 
WN 
128) 
rat 
0) 
N 
© 
iat 


ha 


LIL L/L f | 
LULL LLL. 
HETTTES 


SIMO TEE 
LLL. 


Ge 


Mh. 
a aa yi 


yy) 


Lay 


SL LL 


of 


LLL L110 
LLL LPL 


an SY 
© 
wo 


aa 
A A 


ES eS 
[f/f /[/ 


TED | 


H-OIl€ WNSS3SD 
eulbuy Olb-O] 


LLL ff 1 LLL 


Oo 
aN 
aN 


Oo 
wo 
N 


sHV UG 


(SNIONS 3NO) YSMOdSSYOH 


JOOUDS sfONpDIO {SOY JDADN 
SONMOUOIOY Jo LUdIUuDdsG 


BH, ‘3UNSSSAUd GTIOSINYN 
O€ 62 82 le 92 GS bS SS) ZZ 1ZO7 G! AI ZI 
SS SS Gs es Ss OS ee NO SS | 


SONVWYOSY3d WA31 v3S 





PIS] = wodJ 
2SD3JIUI,9 Ye 


6 | a! 
SO for 7 
: Lo 40} Yt LadaIGng 2 
ge *pIS] WO} osDa299p 


oF YIDS 10) YI ppy'| 


NOILOSYYOD SYNLVYSdW3L 


Pusure 2 


ENGINE PERFORMANCE CURVES 


“1-2 





greatly simplifies the method by which manifold pressure is 
accounted for. Since the altitude curves are graphed as 

a linear function of the density ratio, a simple conver- 
Sion is made from manifold pressure to density ratio and 
cCHismeOnverLed raciom called MAPSIG) is used for calcula— 
tions, as will be shown later. Table I shows the data used 
for the MAPSIG conversion. This is a purely mathematical 
manipulation, used for simplifying the equations by restating 
one variable in terms of another. This converted ratio is 
related to the critical operating altitude and is there- 
Porc sem@icrenu lal com Tone the trie density ratio (SIGMA) ter 
the actual operating altitude, which is determined in the 
program from the flight altitude. Table II shows the data 
used for determining the flight altitude density ratio. 

Once the engine curve parameters have been determined, power 
requirements are calculated for a single engine. 

On the sea level side, brake horsepower (SSLBHP) is 
calculated as a function of the maximum brake horsepower 
(MAXBHP) and maximum manifold ta (MAXMAP) for the 
operating RPM, slope of the RPM line (RSLOPE), and operating 
manifold pressure (MAP) according to the equation: 

SSLBHP = MAXBHP — RSLOPE (MAXMAP - MAP) (1) 

The procedure for determining a value Bow Tomales horsepower 
at altitude is reasonably straightforward. The critical 
altitude (BHPCR) for a given manifold pressure is that 
altitude represented by the intersection of the operating 
RPM and manifold pressure lines. The converted density 


ratio (MAPSIG) discussed earlier Gefines this altitude. 
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TABLE I 


MAPSIG CONVERSION DATA 


MAP MAPSIG MAP MAPSIG MAP MAPSIG 
1h. ou ZO, S70 20 S05 

15. 2580 Pe ~760 wh) De) AOS 

Oy: 0015 ks ~790 Oe 0965 

ie 64,0 pon ~820 Asie -990 

18. OVS 2h O55 

19. 70> Ae 880 

TABLE II - 


ALTITUDE/DENSITY RATIO CONVERSION DATA 


A SIGMA Po ST GMA Niall SIGMA 
Oy ab olelorele 6000. » 83586 ILA. ie) S14 

1000. °97106 OOO, °81064 1300@,) -0/ 133 

2000. 094277 S000. 78601 i} OOOkae e002 

3000. 98512 9000. «76196 M5 COCO O29 23 

1.000. «88808 10000. 073848 

5000. «86167 11000. ~ L S, 





A straight line is drawn connecting this point with the 
sea level brake horsepower value. Then a vertical line is 
drawn upward from the operating density ratio on the abscissa. 
If this line intersects the SSLBHP-—BHPCR line, then that 
intersection defines the brake horsepower required for 
level flight at the operating altitude. If the vertical 
lianem 1es! tComtme right ol the critical, altitude then its 
intersection with the operating RPM line defines the 
required brake horsepower. Figure 3 shows an example of 
engine curve utilization. 

Critacal brake horsepower (BHPOH) is calculated as a 
function of maximum brake horsepower for the operating RPM 
(MABHP), the slope of the operating RPM line (ARSLOP) and 
the converted density ratio (MAPSIG) according to the 


equation: 
BHPCR = MABHP — ARSLOP (1.0 — MAPSIG) (CZ) 


The slope of the line connecting sea level and critical 
brake horsepowers (LSLOPE) is also a function of the con- 


verted density ratio: 
LSLOPE = (BHPCR — SSLBHP) / (1.0 —- MAPSIG) (3) 


The operating altitude is compared with the critical altitude. 
If the operating altitude is the lower one, brake horsepower 


is calculated on the operating RPM line as: 


BHPALT = SSLBHP + LSLOPE (1.0 —- SIGMA) (4) 
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If the operating altitude is above the critical altitude 
then brake horsepower is calculated on the operating RPM 


jine ass: 
BHPALT = MABHP - ARSLOP (1.0 — SIGMA) C5) 


It should be noted at this point that the actual flight 
data are based on thrust horsepower rather than brake horse- 
power. In the classroom the student takes an ideal brake 
horsepower, applies the necessary corrections and then mul- 
tiplies that value by propeller efficiency to get the actual 
thrust horsepower available. This program does exactly the 
opposite. From in-flight RPM and manifold pressure it cal- 
culates the actual thrust horsepower required. This is 
corrected, and then converted to shaft horsepower (it is 
assumed that shaft and brake horsepower are equal) by divid- 
fi? wear Ope ll ermmel Lic. CNC ya. 

Since the engine performance curves are plotted in 
terms of brake horsepower, the nomenclature of the program 
up to this point has been written in terms of brake horse- 
power. ‘However, since the horsepower value just found using 
equation (5) is actually a thrust horsepower, from here on 
it will be referred to as a thrust horsepower, until it is 
converted to shaft horsepower. 

Once the basic thrust horsepower value has been deter- 

mined, several correction factors are applied, the first 
of which is for temperature variation. The engine curves 


are plotted for standard temperature: a correction factor 
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Figure 3 
ENGINE CURVE UTILIZATION 


is supplied with them. First the temperature is changed 
from degrees Fahrenheit (TEMPF) to degrees Centigrade 
(DEGC) and the difference between operating and standard 


temperatures (TDIF) is determined: 
DEGC = 5.0 (TEMPF - 32.0) / 9.0 (6) 


Tonos 150 Smee « (7) 


II 


Then, using the correction factor, the single engine thrust 


horsepower corrected for temperature variation (CBHP) is: 
CBHP = BHPALT (1.0 -0.0016667 (TDIF)) (8) 


The total thrust horsepower (TBHPA) is found by 
multiplying the temperature corrected value by the number 


of engines: 
TBHPA = CBHP (NOENG) (9) 


Total shaft horsepower required at altitude is found 


by dividing thrust horsepower by propeller efficiency: 
TSHPA = TBHPA / PE Gilfe)) 


Propeller efficiency has been arbitrarily chosen as 0.75 
and is entered in the program using a DATA statement for 
case le accesco monoOuldesuoseduentsevelluacion dictate a 
Change in the value. 

The total shaft horsepower is good for one altitude 


and one temperature. To make this information more useful 











it is converted to an equivalent shaft horsepower (SSLPWR} 
which may be used to find horsepower at any altitude for 


the given temperature: 
SSLPWR = TSHPA /SIGMA (11) 


The relationship between velocity and power and their equiv- 
alent counterparts is most easily shown by the following 


sequence of equations: 


2 
‘ea 3 fo uv 
and 
2 
alg 
Gait ~2 FP Y 
also 
te SP gy ee , P= PF, 
SO 
Pe 
dain = Broo V. 
Define equivalent velocity: ue = oy V 
then 
2 
Galt . (ae We 
| gee ae AT Pe — ag 
=e ache Teas SiC y? 
ed fa we Was 
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Define equivalent power: 5 P 


2 oe BSC, 2 (Je Vv)? 
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Once equivalent power has been calculated, the remaining 
step in the power requirements section is to change velocity 


GO equivalent velocity (EQVEL): 
EQVEL = VEL /SIGMA : Gp 


GC. LIFT AND DRAG 

IMeae 1S dLuncllLON Oo, density, mevoe #657, characteristic 
aredmancdeume lai tecocil|1Clengg in Order Go avoid the 
necessity of inputing or calculating a density for the 
flight altitude, the equations for determining lift and 
drag coefficients are based on equivalent velocity, reducing 
density to a constant (sea level value). 


Litt. —ewegtte sa woe C. 
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For the calculations, equivalent velocity is changed from 
knots to feet per second (VEFPS): 

VEFPS = EQVEL (6000.) / 3600. (alg) 
Then the lift coefficient (CLIFT) is found as a function 
of weight (WI) and equivalent velocity: 

CLIFT = 2 (WE) / (.0023769) (S) (VEFPS)© (14) 

The characteristic area, S, is defined by a data state- 

ment at the beginning of the program for ease of access. 
There are two values given for S (NASA TN-—D-—6238, 1971). 
In this program 200 square feet is used, which includes tip 


tanks. Without them the wing area is 179 square feet. 


A 








Drag is a AUnctblon Of densaty, veloagnty, characteristic 
area and the drag coefficient. It is also a function of 
velocity and power: 

Drag = P/V =P, Vas 
=¢ EV, SC, 
SE We SC. Vv" va 
3 & 


or equivalently, Vi 


=2/% 5 C6, V, 
For the calculations, equivalent power is changed from 
shaft horsepower to foot—pounds per second (PEFPPS): 
PEEPP oe mooLPWhmo ene fea) 
igen sehiomaracmeocr frewentp (CORAG\e1s foundsaseas funceien 
of equivalent power and equivalent velocity: 
CDRAG = 2 (PEFPPS) / (.0023769) (VEFPS)? (16) 
Once all computations have been completed, a table of 
velocity, manifold pressure, brake horsepower at altitude, 
equivalent velocity, equivalent shaft horsepower, lift 
coeificient and drag coeificient is printed out. Then two 
calls are made to DRAWP, one of IBM’s plotting subroutines 
used in conjunction with a CALCOMP plotter. The first call 
is for a graph of equivalent power versus equivalent velo- 
Ci Gyir Megeccons is for a graph of lift coefficient versus 


drag coefficient. 


we INTERPOLATION SUBROUTINE 
As stated earlier, the heart of the program is the 
data obtained from engine curves for the Cessna. Originally 


in WoaSsmecer ded tChagminmesce data would te used to TOnmbLat ec 
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equations with RPM as the independent variable. Then, 
given an operating RPM, all necessary information could be 
obtained using the equations, eliminating the need for 
table searches and interpolation. Several subroutines 
PrOmeLeM sc scremaliltessubroutine Package, including Lo@Pia. 
CHBFT, and SPLIN, were utilized in an effort to formulate 
equations which would accurately represent the engine data. 
However, due to problems with convergence, extended pre- 
cision requirements for curve coefficients, and decreasing 
accuracy with increasing curve order, all of these routines 
proved unsatisfactory in producing the required equations. 
It was therefore decided to write a simple subroutine which 
would interpolate in a table of data from the engine curves. 
oubroutine LINT performs linear interpolation or 
extrapolation. Inputs to the subroutine are two columns 
of data and an independent variable for which LINT will 
interpolate the dependent variable. The data column 
corresponding to the independent variable must be in 
monotonically increasing order but may be arbitrarily 
Spaced. Since the engine curves are for a single engine, 
flight with two engines allows operation of each engine at 
values below those shown on the engine performance curves, 
but not above. Subroutine LINT was therefore written to 
extrapolate below the limits of the curves but not above. 
In the event that the independent variable falls above 
the range of data given, an error message so stating this 


1S printed out. 


Zo 





The method used by LINT is Wee Simple. First a check 
is made to see if the variable is bYSllion: the range for inter- 
polation; if this is so, an extrapolation is performed. 
Otherwise the column is searched until bracketing values 
are found, at which time an interpolation is AO - ise 
a bracketing pair Panne be found, the variable is above 


the data range and the error message is printed. 


Be SORTING SUBROUTINE 

To make the program easier to use, it was written so 
that the student’s data may be loaded in any order. It is 
not necessary, for example, to put the data in order of 
increasing velocity. Subroutine SORT does this ordering 
during program execution so that the printed output data 
is ordered in terms of increasing velocity. SORT is 
written to evaluate one column of data, and to put two 
associated columns of data in the same chronological order 
as the one being evaluated. 

The method used by SORT is simple. A logical flag 
is set equal to FALSE. Then sequential pairs of values 
in the first data column are compared. If the second 
value in any one of these pairs is less than the first 
value in the pair, the positions of the two values, and 
of the associated values in the other two columns, are 
Switched and the logical flag is set equal to TRUE. 
When the entire column has been searched, the flag is 
Giteekca nmin i> las weemercsey to TRUE), this indicates 


the possibility that there may be more switches necessary 


ek 





so the flag is reset to FALSE and the procedure is repeated. 
When a complete column search is made without resetting 

the flag, the data are in proper order and control returns 
to the main program. Like the main program, the subroutine 


fowtimiceds co 100 data points per column. 


IIT. PROGRAM USAGK 


Ae DATA LOADING 

The first data to be loaded following execution of 
the program are the engine data, consisting of the slope, 
maximum manifold pressure and maximum brake horsepower 
(again, brake horsepower values are used here because the 
Eneinesdavascomes trom curves which are plotted im terms 
of brake horsepower) for sea level performance, and slope 
and maximum brake horsepower for altitude performance, for 
each of nine sets of RPM lines. These data are shown in 
Table III. Following the permanent data are the students 
in-flight data. 

The order in which in-flight data is loaded varies. 
The program gives the option of constant or variable RPM, 
and the input order depends on which option is selected. 
To simplify the student data loading, input data formats 
for both options are included in the program. These formats 
are similar to those used by James Melsa in his book of 
computer programs for linear control theory. They give 
card by card instructions for loading the data, including 


card number, column numbers, data description, and input 
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RPM 


1800 
1900 
2000 
2100 
2200 
2300 
24.00 
2500 
2625 


RPM 


2100 
2100 
2100 
2100 
2200 
2300 
2300 
FarAele) 
2100 
2100 
2100 
FAG), 
214.50 


TABLE III 


ENGINE PERFORMANCE DATA 


MA BHP 
CHP ) 


5 Ore © 
1630 
1. Sea 
191.0 
2075 
2A oe 
Caw 
249 20 
260.0 


2300 


26 


RSLOPE MAXMAP MAX BHP ARSLOP 
(BHP/INCH) (INCHES HG) (HP) (BHP/INCH) 
OneS 29.75 150.0 173.91 
7016 DOO 164.5 186 .67 
(oem 29.65 ey ae 200.00 
Sxl 7 29.60 3, 0 AMOS) 7/ 
Siew 2 25.55 208.5 PUBS GANS 
92h 29-90 Saat 25h Ih 
10.00 29 4.0 239.0 267 .03 
Ges 5 29 h.0 2148.0 Zoi oe 
Tees 29335 260.0 29h 50 
TABLE LV 
EXAMPLE PROBLEM DATA 

MAP VEL RPM MAP 
(INCHES HG) (KNOTS) (INCHES HG) 
1 850 24.50 20.0 
19.0 100.0 2300 18.0 
O75 11540 2100 18.0 
20.8 30.0 2100 SIC 
2EA@ Wh, 0 2100 TAS. 
AIS AS 160.0 2100 15730 
230 160.0 2100 LEAS 
2240 Wigs 8, 2100 So 
20-i i020 2100 19.0 
I 2 ee 2100 age 
ILS ale OG 40 2300 Ze 
IM 65.0 2300 DSR. 
PS) SO. 161.0 PARC. 


VEL 


(KNOTS) 


ele 
135.0 
Ie 
10626 
8765 
80.0 
S90 
Jer 9, 
Wao) 
Oe 
14h.) 
are 
5S 30 








format. These formats are presented in the program in 
Appendix A. 

Following the input of constant flight parameters, 
such as altitude and weight, velocities and manifold pres- 
sures are read ine If RPM is variable, its values¥are 
ead an; otherwise the constant RPM’ valuwi@ is read iin and 
duplicated for each set of velocity and pressure. This 
duplication is done for use in the sorting subroutine, 
which requires a column of RPM data, and for the printed 
output, which includes a column for RPM. Once all the 
data are read in and sorted, the engine curve performance 
parameters are interpolated for each RPM, if RPM is vari- 
alc. ii aie iS Constanl® the parameters aremevaluated only 
once, saving unnecessary computer time. Following this 
evaluacton, the power, Jiit and drag requirements are 


ee kel ated. 


B. GRAPHICAL OUTPUT 

The power requirements and lift and drag coefficients 
for the aircraft are graphed on the CALCOMP plotter, using 
TBM*s DRAWP subroutine. Each graph is eight inches by 
ten inches, with a title below. The abscissa and ordinate 
are automatically scaled by the subroutine and the scale 
information is printed below the graph. Data points are 
plotted on the graph as triangles. This was done for 
several reasons, If the subroutine is instructed to plot 
a continuous curve, it merely connects each point to the 


next with a straight line. This will produce a very 


Al 





erratic curve in the case of data with any degree of spread, 
and where there are only a few data points to plot. The 
use of triangles allows the student to draw a smooth curve 
through the data. It also allows the data to be loaded in 
any order, rather than forcing the student, for example, to 
arrange his data in increasing order. This last point is 
admittedly academic, since addition of*cthe sorting sub- 
routine eliminates any plotting problems which might have 


Securreagurom unordered datas 


IV. EXAMPLE PROBLEM 


Two actual sets of student data were combined for the 
example problem. RPM was variable. Table IV shows the RPM 
values and associated manifold pressures and velocities. 
‘Temperature varied slightly during the test flights; an 
average value was calculated for use in the program. Air- 
craft and flight parameters are shown in Table V. The 
order of the data deck is shown in Appendix B. 

Results of the program compared very favorably with 
hand calculated data. Table VI shows the tabulated output 
from the program. There was some spread in the data, which 
is easily visible in the sample graphical output shown in 


Appendix C. 
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TABLE V 
BXAMPLE PROBLEM - AIRCRAFT FLIGHT PARAMETERS 


Number of engines operating: 2 
Flight altitude (feet): 5000 

Temperature (degrees F): 76.55 
Aircraft weight (pounds): 4500 


Number of data points: 26 
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Ve. PROJECT EVALUATION 


It is felt that the program has met the objectives 
previously set forth. A simple, easy-to-use program has 
been written which evaluates aircraft in-flight data and 
produces graphs of power versus velocify, and lift coeffi- 
client versus drag coefficient. It yields results which are 
comparable in accuracy to results calculated by hand. 
Additionally, there is a considerable savings in man-hours 
required by the student to evaluate his data, which was the 


original impetus for this thesis. 


VI. AREAS FOR PROGRAM EXPANSION 


There are several possibilities for future expansion 
of the program. Temperature variation, though small, may 
have a noticeable effect on performance figures, and a pro- 
gram option allowing variable temperatures would be bene- 
ficial. The same is true for an option which would allow 
altitude to vary. Test flights are subject to weather, and 
clouds, for example, do not always allow a flight to be 
performed at a single altitude. Another variable factor is 
weight. Fuel accounts for nearly 16 percent of the aircraft 
maximum gross weight, and on a long flight this could 
appreciably alter performance data. An option which uti- 


lizes fuel readings might well be a worthwhile addition. 


Sa 





Another section which might be added to the program ale: 
one which would plot rate of climb curves for the aircraft. 
This could be a simple subroutine which would pernine inputs 
of time to climb through a given ruee Cmanpre Oxen md: 
range of velocities. An easy conversion to feet per minute 
could be made, and the results plotted. 

Other future possibilities includa the area of dynamic 
stability and control information, such as roll rates, which 


are covered in the second half of the flight evaluation 


techniques course sequencee 
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ENGINE DATA FOR USE IN INTERPOLATION 
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APPENDIX B 


SAMPLE INPUT DATA DECK 
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APPENDIX D 
PROGRAM GLOSSARY 


ALT - (input) Altitude (feet) 
ARSLOP - Slope of any altitude performance RPM line 


BHPALT —- Single engine uncorrected thrust horsepower at 
altitude 


BHPCR ~ Critical horsepower at altitude for any given RPM 
and MAP 


Cimewconyersuen tacbor tor changing velocity to fect per 
second 


C2 — Constant for use in determining drag coefficient 

hoe Gcnstanv tor use air determining litt coefficient 

CBHP - BHPALT corrected for temperature variation 

CDRAG — Drag coefficient 

CLIFT - Lift coefficient 

DEGC —- Temperature, degrees Centigrade 

DUMMY -— Dummy variable used for temporary storage 

EQVEL - Sea level equivalent velocity (knots) 

TALT -— Integer version of ALT 

IPTS - Integer version of POINTS 

ITB - Vector of constants for use in DRAWP subroutine 

IW! —- Integer version of WT 

LINT - Linear interpolation subroutine 

LSLOPE - Slope of the line connecting sea level and altitude 
degree BHP (SSLBHP and BHPCR) for any given RPM 


MABHP - Maximum horsepower for any given altitude performance 
RPM line 


MAP -— (input) Manifold pressure (inches Hg) 
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MAPSIG - Performance curve equivalent value of density ratio 
for any given value of MAP 


MAXBHP - Maximum horsepower for any given sea level perfor- 
mance RPM line 


MAXMAP - Maximum manifold pressure for any given sea level 
performance RPM line 


NOENG — (input) Number of engines operating 

ONERPM - (input) Constant/variable RPM flag 

PE — Propeller efficiency 

PEFPPS - Sea level equivalent power (foot pounds per second) 
POINTS —- (input) Number of data points 


RA -— Vector of altitudes used for converting altitude to 
density ratio 


RC - Matrix of engine performance curve data 
RC1-RC6 - Vector equivalents of data contained in RC 


RMAP - Vector of manifold pressures used for MAPSIG con- 
version 


RPM - (input) Engine revolutions per minute 


RS - Vector of density ratios used for converting altitude 
to density ratio 


RSIG - Vector of density ratios used for MAPSIG conversion 
RSLOPE — Slope of any sea level performance RPM line 

RTB1 - Vector of data used in DRAWP subroutine 

RTB2 —- Vector of data used in DRAWP subroutine 

S . Aircraft characteristic area (square feet) 

SIGMA - Density ratio for a given altitude 

SORT - Sorting oe: 

SSIL.BHP —- Sea level power for any given RPM and MAP 

SSTLPWR - Sea level equivalent total shaft horsepower 


TBHPA - Total horsepower at altitude corrected for tCempera— 
ture and number of engines 
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TDIF - Difference between standard temperature and tem- 
perature at flight altitude 


TEMPF - (input) Temperature at altitude (degrees Fahrenheit) 
TITLE1L - Graph title information for DRAWP subroutine 
TITLE2 - Graph title information for DRAWP subroutine 

TSHPA — TBHPA corrected for propeller efficiency 

VEFPS - Sea level equivalent velocity (feet per second) 

VEL - (input) Velocity (knots) 


Wr - (input) Aircraft weight (pounds) 
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